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The compound CaV2 04 contains V"*"^ cations with spin S = 1 and has an orthorhombic structure 
at room temperature containing zigzag chains of V atoms running along the c-axis. We have grown 
single crystals of CaV204 and report crystallography, static magnetization, magnetic susceptibil- 
ity X, ac magnetic susceptibility, heat capacity Cp, and thermal expansion measurements in the 
temperature T range of 1.8-350 K on the single crystals and on polycrystalline samples. An or- 
thorhombic to monoclinic structural distortion and a long-range antiferromagnetic (AF) transition 
were found at sample-dependent temperatures Ts ~ 108-145 K and Tn « 51-76 K, respectively. 
In two annealed single crystals, another transition was found at ~ 200 K. In one of the crystals, 
this transition is mostly due to V2O3 impurity phase that grows coherently in the crystals during 
annealing. However, in the other crystal the origin of this transition at 200 K is unknown. The xiT) 
shows a broad maximum at ~ 300 K associated with short-range AF ordering and the anisotropy of 
X above Tn is small. The anisotropic x{T — > 0) data below Tn show that the (average) easy axis of 
the AF magnetic structure is the 6-ajcis. The Cp{T) data indicate strong short-range AF ordering 
above Tn, consistent with the xiT) data. We fitted our x data by a J1-J2 S — 1 Heisenberg chain 
model, where Ji(J2) is the (next)-nearest-neighbor exchange interaction. We find Ji ~ 230 K, and 
surprisingly, J2/J1 ~ (or J1/J2 ~ 0). The interaction J± between these S = 1 chains leading to 
long-range AF ordering at Tn is estimated to be J±/Ji > 0.04. 

PACS numbers: 75.40.Cx, 75.50.Ee, 75.10.Pq, 81.10.Fq 



I. INTRODUCTION 

Low-dimensional frustrated spin systems have rich 
phase diagrams arising from a complex interplay of ther- 
mal and quantum fluctuations and competing magnetic 
interactions at low temperatures. While spin S = 1/2 an- 
tiferromagnetic (AF) chainsi and odd-leg ladders^ have 
gapless magnetic excitations, S = 1 chains and 5 = 1/2 
even-leg ladders with nearest-neighbor (NN, Ji) inter- 
actions exhibit a finite energy gap between the ground 
state and the lowest excited magnetic states. However, 
numerical calculations have shown that the influence of 
frustrating next-nearest-neighbor (NNN, J2) interactions 
play a significant role and depending on the J2/J1 ra- 
tio, can lead to incommensurate helical spin structures 
which may be gapped or gapless.— i^i^i^i^ Such a system is 
described by the XX Z Hamiltonian^ 

2 

H — ^ Jp ^(S'f S'/^p + SfSf_^p + \Si Sl\p), (1) 
p=i I 

where S; is the spin operator at the lih. site, Jp is 
the AF interaction between the NN (p = 1) and NNN 
{p = 2) spin pairs, and A is the exchange anisotropy. For 
j = J2I J\ > 1/4, the classical AF chain exhibits incom- 



mensurate helical long-range ordering described by the 
wave vector q = arccos[— l/(4j)] and a finite vector chi- 
rality k = S; x S^+i which describes the sense of rotation 
(left or right handed) of the spins along the helix. In 
the large-j, small- A limit of the 5=1 chain, one finds 
a corresponding phase^ where spin correlations decay, as 
required for a one-dimensional system, although only al- 
gebraically, but chirality is still long-range ordered. This 
phase is called the chiral gapless phase and is seen to ex- 
ist for all spin quantum numbers For smaller j, a 
chiral gapped phase is observed in the 5=1 chain,— with 
chiral long-range order and exponentially decaying spin 
correlations. 

The above chiral phases are ground state phases of a 
spin system. In a related prediction. Villain suggested 
three decades ago that a long-range ordered vector chi- 
ral phase can exist above the Neel temperature Tn of a 
quasi-one-dimensional spin chain system showing helical 
magnetic ordering below Tni^ This chiral phase would 
have a transition temperature Tq > Tn that could be 
detected using heat capacity measurements.?^ 

The compound CaV204, containing crystallographic 
V"'"'^ spin 5=1 zigzag chains, has been suggested as 
a model experimental system to study the above chiral 
gapless phascj ^^d^ CaV204 crystallizes in the CaFe204 
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FIG. 1: End-on (top) and inclined (bottom) views along the 
c-axis of the CaV2 04 structure showing the V zigzag chains 
with the V atoms in distorted edge- and corner-sharing octa- 
hedral coordination by oxygen. 

structure at room tcmperaturei^iii^ with the orthorhom- 
bic space group Pnam and with all the atoms in distinct 
WyckofF positions 4(c) (x,?/,l/4) in the unit cell. As 
shown in Figs. [T] and [21 two zigzag chains of distorted 
edge- and corner-sharing VOg octrahedra occur within 
the unit cell and run parallel to the c-axis, with the Ca 
ions situated in tunnels between the chains. Two sets of 
crystallographically inequivalent V atoms occupy the two 
zigzag chains, respectively. The VOg octahedra within a 
zigzag chain share corners with the octahedra in the ad- 
jacent zigzag chain. Within each zigzag chain, in order 
to be consistent with our theoretical modeling later in 
Sec. IIVI of the paper, the nearest neighbors are defined 
to be those on different legs of the zigzag chain where 
the NN V-V distance is 3.07 A. The NNN V-V distance 
(3.01 A) is defined to be along a leg of the zigzag chain. 
The similarity between these two distances in CaV2 04 
suggests that Ji « J2, which would result in geometri- 
cally frustrated AF interactions in this insulating low- 
dimensional systemiiSiii 

Previous studies on polycrystalline samples of CaV2 04 
have offered contrasting views on the nature of the mag- 




FIG. 2: Top panel: the skeletal structure of CaV204 showing 
the zigzag V-V chains. The large spheres represent V atoms, 
the small dark spheres Ca atoms, and the small light spheres 
O atoms. Middle and bottom panels: cross-sections of the b-c 
and a-b planes, with the V and O atoms labeled as described 
in Table HIl The Ca atoms have been omitted for clarity. 



netic ground state. Magnetic neutron diffraction mea- 
surements on CaV204 (Ref. [H) gave clear evidence for 
the presence of long-range antiferromagnetic ordering at 
4.2 K (the temperature dependence was not studied, and 
the Neel temperature was not determined). A doubled 
magnetic unit cell along the b and c directions was found 
with AF propagation vector (0 ^ i) and three coUinear 
AF models with the V ordered moments parallel to the b- 
axis were considered. Interestingly, the ordered moment 
was found to be 1.06(6) /iB/(V atom), where fiB is the 
Bohr magneton. This value is a factor of two smaller 
than the value gSfiB = 2.0 /ib/(V atom) expected for a 
spin S — 1 with spectroscopic splitting factor (g-factor) 
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g = 2. Magnetic susceptibility measurementsiSiii showed 
a broad maximum at ~ 250 K, indicating the onset of 
strong short-range AF ordering in a low-dimensional spin 
system upon cooling. The data also showed a finite value 
at the lowest temperatures, indicating that an energy 
gap for spin excitations did not occur, consistent with 
the neutron diffraction measurements. However, these 
magnetic susceptibility data also showed a bifurcation 
below '--^ 20 K between low field (100 Oe) zero-field- 
cooled and field-cooled measurements that was sugges- 
tive of a spin-giass-like freezing rather than long-range 
AF ordering. ^^V nuclear magnetic resonance (NMR) 
measurements^iii showed a nuclear spin-lattice relax- 
ation rate 1/Ti oc T at low temperatures from 2 K to 
30 K, of unknown origin, but again indicating lack of an 
energy gap for magnetic excitations. The authors^*'' W. 
proposed a chiral gapless ordered phase at low temper- 
atures in accordance with theoretical predictions for a 
S' = 1 frustrated XY or XXZ chain. The chiral phase 
implies a helical spin arrangement in contrast to the 
coUinear spin models proposedi^ in the neutron diffrac- 
tion study. Furthermore, the observation of a ^^V nu- 
clear resonance at the normal ^^V Larmor frequenc y^^i^^ 
at temperatures at and below 4 K is not consistent with 
the long-range antiferromagnetic ordering found at 4 K 
from the neutron diffraction measurements,^'^ since such 
ordering produces a very large static local magnetic field 
of order 20 T at the positions of the V nuclei. 

In order to resolve the above inconsistencies regarding 
the magnetic ground state of CaV204 and to search for 
interesting physics in this system associated with possi- 
ble geometric frustration within the zigzag spin chains, 
we have for the first time (to our knowledge) grown sin- 
gle crystals of this compound, and report herein crystal 
structure, static magnetization and magnetic susceptibil- 
ity x{T), ac magnetic susceptibility Xa.c{T), heat capacity 
Cp(r), and anisotropic linear thermal expansion ai{T) 
(i = X, y, z) measurements over the temperature T range 
1.8 to 350 K on polycrystalline and single crystal samples. 
Our x(2^) and Xac(T) measurements do not show any 
signature of a spin-glass transition around 20 K as pre- 
viously reportediiSiii We instead observe long-range an- 
tiferromagnetic (AF) ordering at sample-dependent Neel 
temperatures Tn ~ 51-76 K. 

We have recently reported elsewhere the results of 
^^O NMR measurements on a polycrystalline sample of 
CaV2 04 and find a clear signature of AF ordering at 
Tn ~ 78 We find no evidence of a ^^V NMR signal 
at the normal Larmor frequency at temperatures between 
4 K and 300 K, in conflict with the above previous ^^V 
NMR studies which did find such a resonance jiSiiii In 
single crystals, at temperatures below 45 K we do find a 
zero-field ^^V NMR signal where the ^^V nuclei resonate 
in the static component of the local magnetic field gen- 
erated by the long-range AF order below Tn ~ 70 Kfi^ 
The ordered moment at 4.2 K in the crystals was es- 
timated from the zero-field ^^V NMR measurements to 
be 1.3(3) /^b/(V atom), somewhat larger than but still 



consistent with the value 1.06(6) ^b/(V atom) from the 
above neutron diffraction measurements An energy 
gap A for antiferromagnetic spin wave excitations was 
found with a value A/fce — 80(20) K in the tempera- 
ture range 4-45 K, where k-g, is Boltzmann's constant. 
This energy gap was proposed to arise from single-ion 
anisotropy associated with the S = \ Y'^^ ion. A model 
for the antiferromagnetic structure at 4 K was formulated 
in which the magnetic structure consists of two substruc- 
tures, each of which exhibits collinear antiferromagnetic 
order, but where the easy axes of the two substructures 
are at an angle of 19(1)° with respect to each other. The 
average easy axis direction is along the 6-axis, consis- 
tent with our magnetic susceptibility measurements to 
be presented here, and with the easy-axis direction pro- 
posed in the earlier neutron diffraction measurements!^ 
Our magnetic neutron diffraction studies of the antifer- 
romagnetic structure of CaV2 04 single crystals are qual- 
itatively consistent with the NMR analyses; these results 
together with high-temperature (T < 1000 K) magnetic 
susceptibility measurements and their analysis are pre- 
sented elsewhere 

We also find that CaV204 exhibits a weak orthorhom- 
bic to monoclinic structural distortion upon cooling be- 
low a sample-dependent temperature Ts — 108-147 K, 
discovered from our neutron and x-ray diffraction mea- 
surements to be reported in detail elsewhereii^ In our two 
annealed single crystals only, anomalies in the heat ca- 
pacity and thermal expansion are also found at Tsi « 
200 K. From high-energy x-ray diffraction measure- 
ments reported here, we find that in one of the crystals 
the anomaly is most likely primarily due to the metal- 
insulator and structural transitions in V2O3 impurity 
phase that grows coherently in the crystal when it is an- 
nealed. In the other crystal, we still find a small anomaly 
in the heat capacity at Tgi but where the transition in 
the V2O3 impurity phase is at much lower temperature. 
Hence we infer that there is an intrinsic transition in 
our two annealed CaV204 crystals at Tgi with an un- 
known origin. We speculate that this transition may be 
the long-sought chiral ordering transition envisioned by 
Villain^ that was mentioned above. 

From our inelastic neutron scattering results to be pub- 
lished elsewhere)^ we know that the magnetic character 
of CaV204 is quasi-one-dimensional as might be inferred 
from the crystal structure. The largest dispersion of the 
magnetic excitations is along the zigzag chains, which is 
along the orthorhombic c-axis direction, with the disper- 
sion along the two perpendicular directions roughly a fac- 
tor of four smaller. This indicates that the exchange in- 
teractions perpendicular to the zigzag chains are roughly 
an order of magnitude smaller than within a chain. We 
therefore analyze the magnetic susceptibility results here 
in terms of theory for the S — 1 J1-J2 linear Heisenberg 
chain, where Ji{J2) is the (next-)nearest-neighbor inter- 
action along the chain. With respect to the interactions, 
this chain is topologically the same as a zigzag chain 
where Ji is the nearest-neighbor interaction between 
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spins in the two different legs of the zigzag chain and J2 is 
the nearest-neighbor interaction between spins within the 
same leg. We utilize exact diagonalization to calculate 
the magnetic susceptibility and magnetic heat capacity 
for chains containing up to 14 spins S — 1, and quantum 
Monte Carlo simulations of the magnetic susceptibility 
and heat capacity for chains of 30 and 60 spins. We ob- 
tain estimates of Ji , J2 / Ji , the temperature-independent 
orbital susceptibility xo: and the zero-temperature spin 
susceptibilities in CaV204 from comparison of the theory 
with the experimental xi'^) data near room temperature. 
Remarkably, we find that one of the two exchange con- 
stants is very small compared to the other near room 
temperature, as opposed to J1/J2 ~ 1 that is expected 
from the crystal structure. Thus, with respect to the 
magnetic interactions, the zigzag crystallographic chain 
compound acts instead like a linear 5 = 1 Haldane spin 
chain compound. In Ref. we propose that partial or- 
bital ordering is responsible for this unexpected result, 
and su gges t a particular orbital ordering configuration. 
In Ref. U^, we also deduce that below Tg ~ 150 K, the 
monoclinic distortion results in a change in the orbital 
ordering that in turn changes the nature of the spin in- 
teractions from those of a Haldane chain to those of a 
S — 1 two-leg spin ladder. Here we also compare the 
theory for the magnetic heat capacity with the results of 
our heat capacity experiments. We estimate the coupling 
J± between these chains that leads to the long-range AF 
order at Tn. 

The remainder of this paper is organized as follows. 
The synthesis and structural studies are presented in 
Secini The magnetization, magnetic susceptibility, heat 
capacity and thermal expansion measurements are pre- 
sented in Sec. IIIII In Sec. IIVI we consider the origin of 
the heat capacity and thermal expansion anomalies at 
Tgi « 200 K. We then analyze the x(^) data in terms of 
the predictions of exact diagonalization calculations and 
quantum Monte Carlo simulations to obtain Ji, J2/J1 
and xo- Using the same Ji and J2/J1 parameters, we 
compare the predicted behavior of the magnetic heat ca- 
pacity with the experimentally observed heat capacity 
data. We also obtain an estimate of the interchain cou- 
pling J_L giving rise to long-range AF order at Tn. A 
summary of our results is given in Sec. |Vl 



II. SYNTHESIS, SINGLE CRYSTAL GROWTH, 
AND CRYSTAL STRUCTURE OF CaV2 04 

A. Synthesis and Crystal Growth 

Polycrystalline CaV204 was synthesized via solid state 
reaction by first mixing V2O3 (99.995%, MV Labs) with 
CaCOa (99.995%, Althaea Chemicals) or CaO obtained 
by calcining the CaCOa at 1100 °C. The chemicals were 
ground inside a He glove-box, pressed and sintered at 
1200 °C for 96 hours in flowing 4.5%H2-He, as well as 
in sealed quartz tubes when using CaO, with interme- 



diate grindings. Phase purity was confirmed by pow- 
der x-ray diffraction (XRD) patterns obtained using a 
Rigaku Geigerflex diffractometer with Cu Ka radiation 
in the angular range 29 = 10-90° accumulated for 5 s 
per 0.02° step. Thermogravimetric analysis (TGA) at 
800 °C using a Per kin Elmer TGA 7 was used to check 
the oxygen content by oxidizing the sample to CaV206- 
A typical oxygen content of CaV2O3.98io.05 was deter- 
mined, consistent with the initial stoichiometric compo- 
sition CaV204. 

CaV204 was found to melt congruently in an Ar arc 
furnace with negligible mass loss by evaporation. There- 
fore crystal growth was attempted by pulling a crystal 
from the melt in a triarc furnace (99.995% Ar) using a 
tungsten seed rod.^^ The triarc furnace was custom made 
for us by Materials Research Furnaces, Inc. Using 15- 
20 g premeltcd buttons of CaV204, pulling rates of 0.2- 
0.5 mm/min were used to grow ingots of about 3-6 mm 
diameter and 3.0-4.7 cm length. The length of the ingot 
was limited by contraction of the molten region as power 
was lowered to control the crystal diameter. Obtaining 
a single grain was difficult because small fluctuations in 
the arcs coupled with high mobility of the CaV204 melt 
easily caused nucleation of new grains. Out of multi- 
ple growth runs, a reasonably large single grain section 
could be cut out of one of the ingots. The as-grown in- 
got from the triarc furnace and a single crystal isolated 
and aligned from it are shown in Figs, [^a) and (b), re- 
spectively. Due to the tendency for multiple nucleations 
in the triarc furnace, an optical floating zone (OFZ) fur- 
nace was subsequently used for crystal growthji^ Growth 
rates and Ar atmosphere flow rates were optimized to 
successfully grow large crystals of 4-5 mm diameter and 
4-5 cm length starting from sintered polycrystalline rods 
with masses of 8-10 g. An as-grown rod from the OFZ 
furnace is shown in Fig. ^c). 

Powder XRD of crushed sections from the triarc grown 
ingots as well as from the OFZ grown crystals showed sin- 
gle phase CaV204. Laue x-ray diffraction patterns of a 
single-grain section confirmed its single-crystalline char- 
acter and the crystal was found to grow approximately 
along its crystallographic c-axis. The crystals were ori- 
ented and cut to obtain faces aligned perpendicular to 
the principal axes. They were measured as grown (only 
for the triarc grown crystals) as well as after annealing 
in flowing 5%H2-He gas at 1200 °C for up to 96 hours. 



B. Powder and Single Crystal X-ray and Neutron 
Diffraction Measurements 

Rietveld refinements of the powder x-ray diffraction 
data obtained at room temperature were carried out 
using the program DBWS9807a..i2, The refined powder 
XRD patterns from a polycrystalline sample and crushed 
pieces of the triarc and optical floating zone grown crys- 
tals are shown in Fig. [Hand the reflnement results are pre- 
sented in Tables mini and [ml The XRD of the powdered 
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FIG. 3: Single crystals of CaV204 grown using (a), (b) a 
triarc furnace (a numbered division on the scale is 1 mm) and 
(c) an optical floating zone furnace (compared with a U.S. 
penny). 



tural phase transition at temperature Tg from the high 
temperature orthorhombic structure to a low tempera- 
ture monoclinic structure was foundii^ This transition is 
reflected in Fig.Oby a sudden change in the (0 3 1) Bragg 
peak intensity which occurs at a temperature Tg ~ 112 K 
for the as-grown triarc crystal, and « 141 K and « 147 K 
for the annealed triarc and OFZ-grown crystals, respec- 
tively. Due to twinning the orthorhombic (0 3 1) reflec- 
tion splits below the structural phase transition into the 
(0 3 1) and (0 3 1) monoclinic reflections. The total inte- 
grated intensity at this position increases at Tg because 
of the increased mosaic which results in a reduction of 
the extinction effect. The peak in the intensity at 105 K 
for the as-grown triarc crystal is an experimental arti- 
fact due to multiple scattering. The lattice parameters 
of the low temperature monoclinic phase differ very little 
from the orthorhombic phase and the monoclinic angle 
a « 89.3° is close to 90°. This result and the smoothly 
varying signatures in the thermodynamic properties sug- 
gest that the structural transition is of second order and 
involves a small distortion of the orthorhombic structure. 
Full details of the neutron and x-ray diffraction structural 
measurements and results will be presented elsewhereJ^ 
A higher temperature anomaly in the temperature de- 
pendence of the lattice parameters of a powder sample 
was observed by x-ray diffraction over a temperature 
range of 175-200 K. This transition with Tsi ~ 200 K 
was also observed in the magnetic susceptibility, thermal 
expansion, and heat capacity measurements of two an- 
nealed single crystals as described in Sec. IIIII below. In 
the next section we investigate whether there is a struc- 
tural aspect to this phase transition. 



High Energy X-ray Diffraction Measurements 
on Annealed CaV2 04 Single Crystals 



annealed single crystal samples showed a trace amount 
(^ 1-2 mol%) of V2O3. It is curious that no trace was 
found of V2O3 impurity phase in the as-grown crystals, 
and that this impurity phase only formed after annealing 
the crystals. From the refinement results, the structural 
parameters remain relatively unchanged between poly- 
crystalline samples prepared by the solid state route and 
both as-grown and annealed single crystals grown from 
the melt. 

Single-crystal neutron diffraction data were collected 
on the four-circle diffractometer E5 at the BERII reac- 
tor of the Helmholtz-Zentrum Berlin, Germany. A py- 
rolytic graphite monochromator was used to select the 
neutron wavelength A = 2.36 A. Second order contam- 
ination was suppressed below 10"'^ of first order by a 
pyrolytic graphite filter. Bragg reflections of CaV2 04 
were measured with a two-dimensional position sensitive 
•^He detector, 90 x 90 mm^ in area. The sample was 
mounted in a closed-cycle refrigerator, where the temper- 
ature was controlled between 290 K and 6 K. A struc- 



In order to unambiguously determine the crystallo- 
graphic structure of CaV204 at various temperatures, 
to characterize structural phase transitions, and to check 
the crystal perfection, high-energy x-ray diffraction mea- 
surements {E = 99.43 keV) using an area detector were 
performed on two annealed single crystals at the Ad- 
vanced Photon Source at Argonne National Laboratory. 
At this high energy, x-rays probe the bulk of a crystal 
rather than just the near-surface region and, by rocking 
the crystal about both the horizontal and vertical axes 
perpendicular to the incident x-ray beam, an extended 
range of a chosen reciprocal plane can be recorded.^*' For 
these measurements, a crystal was mounted on the cold- 
finger of a closed-cycle refrigerator surrounded by the 
heat shield and vacuum containment using Kapton win- 
dows to avoid extraneous reflections associated with Be 
or the aluminum housing. Two orientations of the crys- 
tal, with either the [001] or [100] direction parallel to 
the incident beam, were studied allowing the recording 
of the (hkO) or {Qkl) reciprocal planes. For each data 
set, the horizontal angle, fj,, was scanned over a range of 
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TABLE I: Structure parameters at room temperature for CaV2 04 forming in the CaFe2 04 structure, refined from powder 

XRD data. Space Group: Pnam (#62); Z = 4; Atomic positions: 4(c), {x,y, 1/4); Profile: Pseudo-Voigt. The overall isotropic 
thermal parameter B is defined within the temperature factor of the intensity as g-^ssm e/\ 



Sample 


Synthesis 


a (A) 


6(A) 


c(A) 


B (A^) 


-Rwp(%) 


i?p (%) 


an-2-116 


1200 °C soUd state 


9.2064(1) 


10.6741(1) 


3.0090(1) 


2.33(4) 


11.17 


8.24 


an-2-50 


Triarc as grown 


9.2241(11) 


10.6976(13) 


3.0046(4) 


1.72(5) 


12.05 


9.27 


an-2-50 


Triarc annealed 1200 °C 


9.2054(3) 


10.6748(3) 


3.0042(1) 


1.58(5) 


14.54 


10.95 


an-3-074 


OFZ annealed 1200 °C 


9.2089(2) 


10.6774(3) 


3.0067(1) 


0.75(5) 


17.46 


12.77 



TABLE II: Atomic positions (x, J/, 1/4) for CaV204 obtained by Rietveld refinement of powder XRD data at room temperature 
for four samples. 



Sample Number 
Synthesis 


an-2-116 

Sofid State 
(1200 °C) 
a;, y 


an-2-50-cl 

Triarc 
(as grown) 
X, y 


an-2-50-cl 

Triarc 
(annealed) " 
X, y 


an-3-074 

OFZ 
(annealed)" 
X, y 


Ca 


0.7550(3), 0.6545(2) 


0.7562(4), 0.6536(3) 


0.7542(4), 0.6544(3) 


0.7536(4), 0.6550(3) 


VI 


0.4329(2), 0.6117(1) 


0.4320(3), 0.6120(2) 


0.4336(3), 0.6120(2) 


0.4331(3), 0.6114(2) 


V2 


0.4202(2), 0.1040(1) 


0.4204(3), 0.1041(2) 


0.4200(3), 0.1043(2) 


0.4209(3), 0.1043(2) 


01 


0.2083(6), 0.1615(5) 


0.2128(8), 0.1593(7) 


0.2049(8), 0.1603(8) 


0.2074(9), 0.1635(9) 


02 


0.1176(5), 0.4744(5) 


0.1157(7), 0.4745(8) 


0.1144(7), 0.4756(8) 


0.1181(9), 0.4738(9) 


03 


0.5190(7), 0.7823(5) 


0.5153(11), 0.7812(7) 


0.5166(0), 0.7806(7) 


0.5169(11), 0.7797(8) 


04 


0.4203(6), 0.4270(5) 


0.4207(8), 0.4282(7) 


0.4244(8), 0.4325(7) 


0.4280(9), 0.4251(9) 



"Two-phase sample containing ~ 1% V2O3 as determined from 
Rietveld refinement. 



TABLE III: Bond angles V-O-V and bond lengths V-V for CaV204 at room temperature refined from powder XRD data 
and calculated using Atoms for Windows, version 5.0. The Vi-O and V2-O bond lengths varied from 1.92 A to 2.08 A. The 

accuracy of the bond angles calculated is ±0.1 °. 



Sample Number 


an-2-116 


an-2-50-cl 


an-2-50-cl 


an-3-074 


Synthesis 


Sintered powder 


Triarc 


Triarc 


Optical float zone 




(1200 °C) 


(as grown) 


(annealed) 


(annealed) 


Vi-Oi-Vi (NN)(°) 


93.9 


92.9 


95.0 


93.7 


V1-O4-V1 (NN)(°) 


93.6 


93.0 


94.4 


96.6 


Vi-Vi (NN)(A) 


3.009 


3.005 


3.004 


3.004 


V1-O4-V1 (NNN)(°) 


99.3 


99.9 


101.8 


100.3 


Vi-Vi (NNN)(A) 


3.078 


3.094 


3.077 


3.071 


V2-O2-V2 (NN)(°) 


93.1 


93.0 


93.7 


92.5 


V2-O3-V2 (NN)(°) 


96.8 


95.6 


95.6 


97.5 


V2-V2 (NN)(A) 


3.009 


3.005 


3.004 


3.004 


V2-O3-V2 (NNN)(°) 


97.3 


98.0 


98.3 


97.1 


V2-V2 (NNN)(A) 


3.058 


3.062 


3.062 


3.055 


V1-O1-V2 {°Y 


121.7 


122.9 


121.8 


121.0 


V1-V2 (A)" 


3.583 


3.581 


3.582 


3.589 


V1-03-V2 (°)'' 


131.6 


132.2 


132.2 


132.2 


V1-V2 (A)' 


3.647 


3.652 


3.652 


3.643 



"Interchain angles and distances parallel to a-axis 
'interchain angles and distances parallel to 6-axis 
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FIG. 4: Rietveld refinement of room temperature powder XRD data of CaV204 showing lobs (+) Icaic {■), difference (-), and 
peak positions (|) for (a) a polycrystalline sample, (b) an as-grown triarc crystal, (c) an annealed triarc crystal, and (d) an 
annealed optical floating zone (OFZ) crystal. The annealed single crystal samples contain small XRD peaks from 1-2 mol% 
of V2O3, shown as the lower sets of peak position markers. 
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FIG. 5: Temperature variation of intensity of the (0 3 1) struc- 
tural Bragg peak across the orthorhombic to monoclinic struc- 
tural transition temperature (Ts) in single crystal samples of 
CaV204 measured by neutron diffraction. The (0 3 1) peak 
is present in both the orthorhombic and monoclinic phases. 



±2.4 deg for each value of the vertical angle, r], between 
±2.4 deg with a step size of 0.2 deg. The total exposure 
time for each frame was 338 sec. The x-ray diffraction 
patterns were recorded with different intensities of the 
incident beam that were selected by attenuation to in- 
crease the dynamic range to a total of 10^ counts. A 
beam size of 0.3 x 0.3 mm^ was chosen to optimize the 
intensity/resolution condition and to allow probing dif- 
ferent sections of the crystal by stepwise translations of 
the crystal in directions perpendicular to the incident 



beam. 



1. Annealed Triarc- Grown Crystal 

Figure [S] shows the (hkO) diffraction plane of the an- 
nealed triarc sample an-2-50-cl at 205 K. The recipro- 
cal space image reveals well-defined diffraction spots that 
correspond to the "primary" CaV204 lattice, as well as 
spots that can be indexed to an impurity phase inclusion 
of V2O3 coherently oriented with respect to the CaV204. 
No additional reflections were observed. Indeed, we find 
two coherent twins of V2O3 related by an inversion across 
a mirror plane of the CaV204 lattice as depicted by the 
black and white rectangles in the left panel of Fig. [6l By 
comparing the integrated intensities of reflections from 
the two phases, we estimate that V2O3 comprises a vol- 
ume fraction of approximately 1-2 percent of the sample. 
This is in excellent agreement with the result of the x- 
ray diffraction analysis of the polycrystalline sample pre- 
pared from the same annealed crystal that was described 
above. The volume fraction of V2O3 varies only slightly 
in different parts of the crystal probed by scanning the 
x-ray beam over the crystal. This indicates that the in- 
clusions of V2O3 are approximately equally distributed 
over the volume of the crystal. 

Upon lowering the temperature of the crystal to 
165 K, below Tgi ~ 200 K, we observe changes in 
the V2O3 structure consistent with the known first- 
order rhombohedral-to-monoclinic structural transition 
at 170 K (measured on heating). ^^'^^ In particular, the 
upper left corners of both panels of Fig. [5] show enlarged 
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FIG. 6: High-energy x-ray difFraction patterns of the annealed triarc CaV204 crystal (an-2-50-cl), oriented with the (ftfcO) 
plane coincident with the scattering plane at T = 205 K (left panel) and 165 K (right panel). The white circles in the lower 
center of each pattern depict the excluded areas around the primary x-ray beam direction. Several peaks corresponding to 
the main phase, CaV204, as well as the coherently oriented second phase, V2O3, are labeled by indices {hkV) and (hkl)* , 
respectively. For V2O3, the hexagonal Miller indices for the rhombohedral lattice are used. The insets of both panels display 
enlarged regions of the diffraction pattern to highlight the rhombohedral-to-monoclinic transition for V2O3. 



views of the region near the (380) reflection from CaV204 
and the (220) reflection (in hexagonal notation) for the 
rhombohedral lattice of V2O3. Below Tgi the (220) re- 
flection splits into three reflections in the monoclinic 
phase. The temperature dependence of this transition 
is displayed in Fig. [T] Here, we note that there is a fi- 
nite range of coexistence between the rhombohedral and 
monoclinic phases of V2O3 (see the inset to Fig. [7]) and 
the transition itself has a hysteresis of roughly 5-10 K. 

Several points regarding Figs. [S] and [7] are relevant to 
our interpretation of the specific heat and thermal expan- 
sion measurements of the annealed triarc crystal (an-2- 
50-cl) to be presented below in Figs. [T51and[TBl respec- 
tively. First, we note that over the temperature range 
encompassing the features at Tsi ^ 200 K, there is no 
apparent change in the diffraction pattern of CaV2 04. 
These anomalies are instead strongly correlated with the 
rhombohedral-to-monoclinic transition in V2O3. We fur- 
ther note that the temperature for this latter transition is 
somewhat higher than the accepted value of « 170 K (de- 
termined on warming) found in the literature. ^"^'^^'^'^ This 
difference is, perhaps, due to the fact that the V2O3 and 
CaV204 lattices are coupled, as evidenced by the coher- 
ent orientation relationship between them, so that strains 
at the phase boundaries come into play and can raise the 
transition temperature^^ In addition, it is report e d^^i^^ 
that deviations of the stoichiometry from V2O3 can affect 
the transition temperature significantly. 

We now turn our attention to changes in the diffrac- 
tion pattern of CaV204 associated with anomalies in the 
heat capacity and thermal expansion measurements at 
temperatures Ts ~ 150 K in Figs. [15] and [16] below, re- 
spectively. The annealed triarc crystal (an-2-50-cl) was 
reoriented so that the CaV204 (Ofc^) reciprocal lattice 
plane was set perpendicular to the incident beam. Fig- 



ure [5] shows the diffraction patterns obtained at 205 K 
(above Ts) and 13 K (well below Ts). The strong re- 
flections in Fig. [5] are associated with the main CaV204 
lattice while the weaker diffraction peaks are, again, asso- 
ciated with the coherently oriented V2O3 second phase. 
At low temperatures, we observe a splitting of the main 
reflections that is the signature of the orthorhonibic-to- 
monoclinic transition at Ts for the CaV204 lattice. The 
splitting of reflections associated with the transition at 
Tsi for V2O3, in this reciprocal lattice plane, is not read- 
ily observed. 

For the low-temperature monoclinic phase of CaV204 
two possible space groups have been considered>i^ The 
space groups P 2i/n 1 1 and P n 1 1 can be separated 
by testing the occurrence or absence of (OfcO) reflections 
with k odd, respectively. The systematic absence of such 
reflections was proven by recording (hkO) planes with 
varying conditions to evaluate the sporadic occurrence 
of these reflections by Rengers or multiple scattering. 
The space group P 2i/n 1 1 is confirmed for the low- 
temperature phase of the studied CaV204 crystal. No 
changes in the diffraction pattern were observed related 
to the onset of antiferromagnetic order in CaV204 below 
Tn = 69 K. 

The details of the orthorhombic-to-monoclinic transi- 
tion at Ts for CaV204 are shown in Fig. [9] where we plot 
the monoclinic distortion angle as a function of temper- 
ature. The monoclinic angle was determined from the 
splitting of the peaks along the b-direction through the 
position of the (042) reflection. Below Ts = 138(2) K, the 
monoclinic angle evolves continuously, consistent with a 
second order transition, and saturates at approximately 
89.2 deg at low temperatures. 
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FIG. 7: Integrated intensity of selected reflections in the high- 
energy x-ray diffraction pattern related to the rhombohedral 
(top panel) and monoclinic (bottom panel) phases of V2O3 as 
a function of temperature for the annealed triarc crystal (an- 
2-50-cl). The direction of the temperature change is indicated 
by arrows. The inset displays the pattern observed at 180 K 
(with increasing temperature) for the (220) reflection from 
V2O3 in the coexistence temperature range. This region is 
the same as that displayed in the insets of Fig. [6] 



13K 



FIG. 8: High-energy x-ray diffraction patterns of the (Okl) 
reciprocal lattice plane of CaV2 04 from the annealed triarc 
crystal (an-2-50-cl) at 205 K (top panel) and 13 K (bottom 
panel). The white circles in the center of the patterns depict 
the excluded areas around the primary x-ray beam direction. 
Most of the reflections related to CaV2 04 show intensities 
above 10^ counts (see intensity scale). The (020) and (002) 
reflections of CaV2 04 are marked in the top panel. The area 
bounded by the black rectangle in the top panel depicts the 
region close to the orthorhombic (042) reflection of CaV204, 
analyzed in Fig. |9l 



2. Annealed OFZ-Grown Crystal 

The annealed optical floating zone crystal (an-3-074 
OFZ) shows a diffraction pattern similar to that of the 
annealed triarc-grown crystal (an-2-50-cl) in measure- 
ments of (hkO) planes at room temperature. The ob- 
served V2O3 inclusions are again coherently oriented with 
respect to the CaV204 lattice. The intensities of the 
diffraction peaks related to V2O3 are similar to those 
in the annealed triarc crystal (an-2-50-cl) and also vary 
only slightly upon scanning different spots of the crystal 
which indicates a homogeneous distribution of the V2O3 
inclusions with a similar volume fraction. However, the 



temperature dependence of the diffraction pattern is dif- 
ferent for the two crystals. Measurements taken on cool- 
ing show that in the annealed floating-zone crystal (an-3- 
074 OFZ), the shape and position of the peaks originat- 
ing from V2O3 are stable from room temperature down 
to 130 K where the onset of the structural transition oc- 
curs. Around 120 K strong changes are observed similar 
to the observations around 180 K in the annealed triarc 
crystal (an-2-50-cl). Below 110 K, the transition to the 
low-temperature monoclinic structure of V2O3 is com- 
plete. Therefore, the temperature for the rhombohedral- 
to-monoclinic transition is reduced by ^ 60 K compared 
to the corresponding temperature in the annealed triarc 
crystal (an-2-50-cl). 
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FIG. 9: Temperature dependence of the monoclinic angle in 
the low-temperature phase of CaV2 04, extracted from diffrac- 
tion data such as shown in the inset for temperatures of 140 K 
and 13 K. Such diffraction line scans were extracted from 
high-energy x-ray diffraction patterns by summing up the in- 
tensity perpendicular to the b direction for the sector marked 
by the rectangle in the top panel of Fig. [S] 



III. MAGNETIZATION, MAGNETIC 
SUSCEPTIBILITY, HEAT CAPACITY AND 
THERMAL EXPANSION MEASUREMENTS 



TABLE IV: Antiferromagnetic ordering (Neel) temperature 
(Tn), high temperature orthorhombic to low-temperature 
monoclinic structural transition temperature (Ts), and the 
transition temperature at ~ 200 K (Tsi) observed by static 
magnetic susceptibility x (peak of d(xr)/dT'), heat capac- 
ity Cp (peak of ACp), thermal expansion a [peak of a{T), 
except for Tsi where the onset of a slope change is used], 
powder synchrotron x-ray diffraction (XRD), single crystal 
neutron diffraction (ND), and single crystal high-energy x-ray 
diffraction (HEXRD) measurements for poly crystalline (pow- 
der) and single crystal CaV2 04 samples. The single crystals 
were grown using either a triarc furnace or an optical floating 
zone (OFZ) furnace. 



Sample Synthesis 


Method Tn (K) Ts (K) Tsi (K) 


an-2-116 powder 


X 


76 


147 


a 


1200 °C 




75 


144 


a 




XRD 




150 


200 


an- 2-50 Triarc crystal 


X 


51 


108 


a 


as- grown 


Cp 


51 


108 


h 




ND 


53 


112 


h 


an- 2-50 Triarc crystal 


X 


68 


133 


195 


annealed 1200 °C 


c-p 


68 


133 


193 




a 


68 


136 


198 




ND 


69 


141 


h 




HEXRD 




138(2) 


192(7)' 


an-3-074 OFZ crystal 


X 


69 


136 


192 


annealed 1200 °C 




71 


132 


191 




ND 


69 


147 


h 



In the following, we describe our results of magneti- 
zation, magnetic susceptibility, heat capacity, and ther- 
mal expansion measurements of both polycrystalline and 
single crystal samples. These and additional measure- 
ments consistently identify temperatures at which the 
antiferromagnetic transition (Tn), the orthorhombic-to- 
monoclinic structural transition (Tg) and the transition 
at ^ 200 K (Tsi) occur in these samples. In Table ITVl we 
summarize these transition temperatures for the different 
samples obtained using the various measurements. 

A. Magnetization and Magnetic Susceptibility 
Measurements 

The static magnetic susceptibility versus temperature 
x{T) = M{T)/H of a polycrystalline sample as well as 
of the oriented crystals was measured using a Quantum 
Design MPMS SQUID magnetometer in a 1 T field from 
1.8 K to 350 K, where M is the magnetization of the sam- 
ple and H is the magnitude of the applied magnetic field. 
In addition, low field (100 Oe) zero- field-cooled and field- 
cooled (zfc, fc) measurements of M{T) at fixed H were 
carried out from 1.8 K to 100 K. A Quantum Design 
MPMS ac SQUID magnetometer was used to measure 
the ac susceptibility Xac(2^) of the annealed triarc grown 
crystal from 5 to 100 K in an ac field _ffac = 1 Oe and 



"No clear signature of a transition corresponding to Tgi was ob- 
served. 

''Not measured for this sample. 

■^The error bar reflects the hysteresis on warming and cooling. 



frequency 10 Hz. The powder was contained in polycar- 
bonate capsules mounted in clear plastic straws. Each 
crystal was glued to a small piece of clear plastic trans- 
parency sheet with GE 7031 varnish or Duco cement, 
which was then aligned inside the plastic straws with the 
a, 6 or c axis direction parallel to the external magnetic 
field. M{H) isotherms were measured in fields up to 
H = 5.5 T at various temperatures. 

The x{T) in = 1 T is plotted in Fig. [TO] for a 
CaV2 04 polycrystalline sample and for aligned single 
crystals grown using a triarc furnace and using an opti- 
cal floating zone (OFZ) furnace. The broad maximum in 
x(T) around 300 K is characteristic of a low-dimensional 
spin system with dominant antiferromagnetic exchange 
interactions with magnitude of order 300 K. For the sin- 
gle crystal samples, clear evidence is seen for long-range 
antiferromagnetic ordering at Neel temperatures Tn = 51 
to 69 K, depending on the sample. The easy axis of the 
antiferromagnetic ordering (with the lowest susceptibil- 
ity as T — > 0) is seen to be the 6-axis, perpendicular 
to the zigzag V chains. At temperatures above Tn, the 
susceptibility of the crystals is nearly isotropic, but with 
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FIG. 10: Magnetic susceptibility x versus temperature in a field of 1 T of CaV2 04 (a) polycrystalline sample, (b) as-grown 
triarc-grown single crystal, (c) annealed triarc-grown single crystal, and (d) annealed OFZ-grown single crystal. The axes {a, b, 
or c) along which the measurements were carried out are as indicated. The insets show d{xT)/dT versus T to highlight the 
transition temperatures. The oscillatory behavior of d(xr)/dr at the higher temperatures, most pronounced in the inset in 
(b), is an artifact generated by the SQUID magnetometer. 
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FIG. 11: Zero-field-cooled and field-cooled susceptibility of CaV204 in a field of 100 Oe measured on (a) polycrystalline powder, 
(b) unannealed triarc grown single crystal (c) annealed triarc grown single crystal, and (d) annealed OFZ grown single crystal. 
Part (c) also shows the ac-susceptibility Xac(r) along the easy 6-axis of the annealed crystal measured in a field ffac = 1 Oe at 
a frequency of 10 Hz. The antiferromagnetic transition temperatures Tn are marked as shown. 



small anisotropies which typically showed Xb > Xa > Xc- 
However, occasionally variations of ±5% in the absolute 
value of Xi{T) were observed between different runs for 
the same crystal axis i that we attribute to sample size 
and positioning effects (radial off-centering) in the second 
order gradiometer coils of the Quantum Design MPMS 
SQUID magnetometer 



The ordering temperatures observed are marked by 
vertical arrows in Fig. [TO] and are highlighted in the plots 
of d{xT)/dT versus T shown in the insets. The various 
transition temperatures are summarized in Table IIVI As 
is typical for a low-dimensional antiferromagnetic system, 
the polycrystalline sample shows only a very weak cusp 
at Tn = 76 K due to averaging over the three principal 
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FIG. 12: Magnetization M versus applied magnetic field H 
isotherms at a temperature of 1.8 K of a polycrystalline sam- 
ple and of an annealed triarc grown single crystal of CaV2 04. 



axis directions, but it is still well-defined as observed in 
the d{xT)/dT vs T plot shown in the inset of Fig. [TUfa). 
In a related study, ^^O NMR measurements on a poly- 
crystalline sample of ^^0-enriched CaV204 gave a clear 
signature of antiferromagnetic ordering at 78 Kfi^ In con- 
trast to the polycrystalline sample, the as-grown crystal 
in Fig. [TUTb) shows a clear and distinct antiferromag- 
netic ordering temperature but with a much lower value 
Tn ~ 51 K. After annealing the crystals. Figs, \10\ c) 
and fTUTd) show that Tn increases to « 69 K, closer to 
that observed in the polycrystalline sample. However, 
the powder average of the annealed single crystal suscep- 
tibility below Tn does not match the susceptibility of the 
polycrystalline sample. The reason for this disagreement 
is unclear at this time. In any case the slow upturn in the 
susceptibility of the powder sample below 40 K in Fig.fTOl 
is evidently intrinsic, due to the powder average of the 
anisotropic susceptibilities, and is not due to magnetic 
impurities. 

The zero- field-cooled (zfc) and field-cooled (fc) x(T) 
measured in a field of 100 Oe for polycrystalline and sin- 
gle crystal samples of CaV204 are plotted in Figs. fTTT aV 
(d). Also shown in Fig. [TlTc) is the real part of the ac 
susceptibility Xac(^) along the easy 6-axis direction of 
the annealed triarc crystal measured in an ac field of 
amplitude 1 Oe at a frequency of 10 Hz. A small irre- 
versibility is observed in Fig. [TT] in all samples between 
the zfc and fc susceptibilities below 30 K. However, 
the Xac(^) measurement in Fig. Illf c) does not show any 
peak in that temperature region, ruling out spin glass- 
like spin freezing which was suggested to occur in powder 
samples from earlier reportsJ^iii The slight irreversibil- 
ity observed may be associated with antiferromagnetic 
domain wall effects. 

In Fig.[T2]we show isothermal M{H) measurements up 
to H = 5.5 T measured at 1.8 K for the polycrystalline 
sample and for the annealed triarc-grown single crystal. 
The behavior is representative of all samples measured. 
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FIG. 13; Heat Capacity Cp versus temperature T in zero 
magnetic field of one polycrystalline sample and three single 
crystal samples of CaV2 04 and of a polycrystalline sample of 
isostructural nonmagnetic CaSc2 04. On this scale, the data 
for the four CaV204 samples are hardly distinguishable. 



We find that M is proportional to H at fields up to at 
least ~ 2 T, indicating the absence of any significant fer- 
romagnetic impurities and the absence of a ferromagnetic 
component to the ordered magnetic structure. 



B. Heat Capacity Measurements 

The heat capacity Cp versus temperature T of a sin- 
tered polycrystalline pellet of CaV204 as well as of 
crystals (as-grown and annealed) was measured using a 
Quantum Design PPMS system at T = 1.8 to 200-300 K 
in zero applied magnetic field. The Cp(T) was also mea- 
sured of a polycrystalline sintered pellet of isostructural 
(at room temperature) nonmagnetic CaSc204 whose lat- 
tice parameters and formula weight are very similar to 
those of CaV2 04.^'* The CaSc204 sample was synthe- 
sized from SC2O3 (99.99%, Alfa) and CaCOg (99.995%, 
Aithaca) by reacting a stoichiometric mixture in air at 
1000 °C for 24 hr and then at 1200 °C for 96 hr with in- 
termediate grindings, and checked for phase purity using 
powder XRD. 

In Fig. [13] we plot the measured Cp{T) of four CaV204 
samples and of isostructural nonmagnetic CaSc204. The 
difference ACp versus T for the four CaV204 sam- 
ples is plotted in Fig. [TD(a). Here ACp is the differ- 
ence between the heat capacity of CaV204 and that 
of CaSc204, but where the temperature axis of Cp 
for CaSc204 was multiplied by a scaling factor to 
take account of the difference in the formula weights 
of CaV204 and CaSc204. This factor is given by 
[MM(CaSc2O4)/MM(CaV2O4)]^/2^0.9705 where Mm is 
the molar mass of the respective compound. If the lat- 
tice heat capacity of CaV204 and the (renormalized) heat 
capacity of CaSc204 had been the same, the difference 
ACp(r) would presumably have been the magnetic heat 
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FIG. 14; (a) ACp versus temperature T for four CaV204 sam- 
ples. Here ACp is the difference between the heat capacity 
of CaV2 04 and that of CaSc2 04, but where the temperature 
axis of Cp for CaSc2 04 was multiplied by 0.9705 to take ac- 
count of the difference in the formula weights of CaV2 04 and 
CaSc204. (b) Entropy AS(r) associated with the ACp(r) 
data in (a), obtained by integrating ACp/T in (a) versus T. 
The AS is normalized by the entropy 2R ln(3) of two moles 
of disordered spins S = 1, where R is the molar gas constant. 
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FIG. 15: ACp versus temperature T for four CaV204 sam- 
ples. The data are the same as in Fig. I14f a) except for ver- 
tical offsets to separate the data sets. The symbol Tn la- 
bels the long-range antiferromagnetic transition and Ts la- 
bels the high temperature orthorhombic to low temperature 
monoclinic structural transition. An additional transition at 
Tsi ~ 200 K is seen in the two annealed single crystals. The 
small anomaly for the AFZ crystal and for the triarc crystal 
that rides on top of the broader peak appear to be intrinsic 
to CaV2 04. The broad peak in the latter crystal appears to 
be due to coherently grown V2O3 impurity phase that grows 
during annealing (see text). 



capacity of CaV204. However, due to the structural tran- 
sition at Ts and the transition at Tsij ^C-p{T) contains 
a lattice contribution as well. The lattice contribution 
to ACp is expected to be minimal at low temperatures, 
where only the long wavelength acoustic phonon modes 
are excited, and possibly also above Tsi where the two 
compounds are known to be isostructural. 

Figure [T4ja) shows that the ACp (T) data for the four 
CaV204 samples are similar except for the different sizes 
and temperatures of the anomalies associated with three 
transitions. In order to more clearly illustrate the dif- 
ferences between samples, Fig. [15] shows the same data 
for each sample but vertically displaced from each other 
to avoid overlap. The magnetic ordering transition at 
Tn as well as the ordering temperatures Ts and Tsi are 
clearly evident from the ACp(T) data in Fig. [15] The or- 
dering temperatures observed are summarized above in 
Table HV] 

The entropy versus temperature associated with the 
ACp(T) data of each sample in Fig. [THa) is shown in 



Fig. [Hb), obtained from AS'(T) = [ACp{T) /T]dT . 
In Fig. [TWb). AS* is normalized by the entropy 
2Rln(2S' + 1) = 2Rln(3) for two moles of fully disordered 
spins = 1, where R is the molar gas constant. As noted 
above, at least at low temperatures, we associate AS{T) 
with the magnetic entropy of the system. At the anti- 
ferromagnetic ordering temperature Tn, the normalized 
value of AS'(TN)/2Rln(3) « 6-8% is very small and is 
about the same for all samples. This small value indi- 
cates that short-range antiferromagnetic ordering is very 
strong above Tn and the data in Fig. [HT b) indicate that 
the maximum spin entropy of the system is not attained 
even at room temperature. This is qualitatively consis- 
tent with our estimate Ji ~ 230 K obtained below in 
Sec. IIVI by comparison of our x(^) data with calcula- 
tions of x(T). 

The small observed magnetic entropy at Tn is consis- 
tent with the values of the heat capacity discontinuities 
ACaf at Tn in Fig.]!^] as follows. In mean field theory, 
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TABLE V: Relative length changes along the orthorhombic 
a-, b- and c-ajds directions and of the volume V that are 
associated with the three transitions at Tn « 70 K, Ts « 
140 K, and Tsi ~ 195 K in annealed triarc CaV204 single 
crystal an-2-50-cl. The temperature T range over which the 
changes were measured are as indicated. 





Tn 


Ts 


Tsi 


r range (K) 


52-75 


120-153 


174-198 


Aa/a (10^^) 


1.2 


4.9 


-7.5 


Ab/b (10^5) 


-4.3 


-0.9 


-0.1 


Ac/c (10-^) 


1.1 


19.6 


-6.3 


AV/V (10-^) 


-2.0 


23.6 


-13.9 



for a system containing A'' spins S the discontinuity in 
the magnetic heat capacity at the ordering temperature 
for either ferromagnetic or antiferromagnetic ordering is 
predicted to be^^ 



5 (25+1)^-1 



(2) 



where N is the number of spins and /cb is Boltzmann's 
constant. Using S = 1 relevant to V"'"'^ and iVfee = 2R, 
where R is the molar gas constant, one obtains ACaf = 
4R= 33.3 J/mol K, where a "mol" refers to a mole of 
CaV204 formula units. From Fig. [151 the experimental 
ACaf is about 0.5 to 2 J/mol K, which is only 1.5-6% of 
the mean field value. This small jump in Cp(Tn) is con- 
sistent with the above small value of 5'(Tn). When short 
range magnetic ordering removes most of the magnetic 
entropy of a system at high temperatures, then thermal 
effects associated with three-dimensional magnetic order- 
ing of the system at low temperatures will necessarily be 
much smaller than otherwise expected. 



C. Thermal Expansion Measurements 

The thermal expansion of the annealed triarc grown 
CaV204 crystal an-2-50-cl was measured versus temper- 
ature using capacitance dilatometry^ from 1.8 to 300 K 
along the three orthorhombic axes a, b and c. The crystal 
is the same annealed triarc crystal measured by mag- 
netic susceptibility and heat capacity in Figs. [10] and 
[TTT c') and in Figs. [TMT51 respectively. In Fig. [TH] the 
linear coefficients of thermal expansion are plotted ver- 
sus temperature (left-hand scale) , along with the volume 
thermal expansion coefficient (right-hand scale) . At high 
temperatures T ~ 300 K the a values tend to become 
temperature independent. Below 200 K, the ordering 
transitions observed above in the magnetic susceptibility 
and heat capacity are reflected in distinct anomalies in 
the thermal expansion coefficients at the corresponding 
temperatures. The ordering temperatures observed are 
summarized above in Table HVl 
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FIG. 16: Coefficients of linear thermal expansion Qa,6,c of an 
annealed single crystal of CaV2 04 measured along the or- 
thorhombic a, b, and c-axes versus temperature (left-hand 
scale) , and the volume expansion coefficient aa + ab + etc ver- 
sus temperature (right-hand scale). 



The normalized length changes along the orthorhom- 
bic a, b and c axis directions and the normalized change 
in the volume V associated with the three transitions at 
Tn, Ts and Tsi are listed in Table [V] These changes 
were calculated by determining the areas under the re- 
spective peaks in the thermal expansion coefficients in 
Fig. [TBI a-iid then subtracting the estimated respective 
background changes over the same temperature intervals. 



IV. ANALYSIS OF EXPERIMENTAL DATA 

A. Origin of the Transition at Tsi ~ 200 K in 
Annealed CaV2 04 Single Crystals 

From Table [y] above, the relative volume change of the 
annealed triarc crystal an-2-50-cl on heating through Tsi 
from the thermal expansion data is AV/V « — 1.4x 10~^. 
This value is about 1% of the value for pure V2O3 at its 
transitioni^ The height of the heat capacity anomaly 
above "background" in Fig. [15] for this crystal is ~ 
2.5 J/mol K, which is about 0.8% of the value^^ at the 
structural transition for pure V2O3. These estimates are 
both consistent with our estimates from x-ray diffraction 
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data in Sec.|TT]of a 1-2 percent volume fraction of V2O3 in 
this crystal. Our data therefore indicate that for the an- 
nealed triarc-grown crystal (an-2-50-cl), the anomalous 
features found above in the heat capacity and thermal 
expansion data at Tgi arise mainly from this transition 
in the V2O3 impurity phase. 

Furthermore, the temperature dependences of the lin- 
ear thermal expansion coefficients at the transition Tgi w 
200 K in Fig. [16] are significantly different than near the 
transitions Tn ~ 70 K and Ts « 140 K. There appears 
to be a discontinuity in the slopes of (T) as the transi- 
tion Tsi is approached from above, whereas a continuous 
change in the slopes occurs as Ts and Tn are approached 
from above. The reason for this difference is evidently 
that the former transition is mainly due to the first or- 
der structural transition in the V2O3 coherently-grown 
impurity phase in this annealed crystal as investigated 
previously in Sec. Ill C H whereas the latter two transi- 
tions are second order. 

However, we also found in SecHIlthat for the annealed 
floating-zone crystal (an-3-074 OFZ), the structural tran- 
sition of the V2O3 impurity phase was reduced by ~ 60 K 
from that of the V2O3 impurity phase in the annealed 
triarc grown crystal, and hence cannot be responsible 
for heat capacity anomaly at Tsi ~ 200 K for the an- 
nealed float-zone crystal in Fig. [151 Indeed, the relatively 
small heat capacity anomaly at Tsi in Fig. [T5| for the 
float-zone crystal appears to also be present at the same 
temperature for the annealed triarc-grown crystal, but 
rides on top of a broader anomaly that is evidently due 
to the structural transition of the V2O3 impurity phase 
in that crystal. Furthermore, the double peak structure 
in the heat capacity for the annealed triarc crystal at 
Tg « 140 K evidently arises due to the overlap of the on- 
sets of the structural transitions in V2O3 and CaV204. 

An issue of interest is the cause(s) of the variabil- 
ity in the structural rhombohedral-to-monoclinic tran- 
sition temperature of the coherently grown V2O3 impu- 
rity phase in our annealed CaV2 04 crystals. Due to the 
first order nature of the transition, the transition is hys- 
teretic. The transition temperature of bulk stoichiomet- 
ric V2O3 has been reported to be at « 170 K on heating 
and « 150 K on cooling i^^i^i^^i^ The transition temper- 
ature decreases rapidly under pressurci^ A pressure of 
only 9 kbar lowers the transition temperature by 60 K, 
and the transition is completely suppressed at a pressure 
of w 20 kbar>2i The transition temperature is also rapidly 
suppressed if the sample contains V vacancies; a crystal 
of composition V1.985O3 showed a transition temperature 
of w 50 K.'^^ Furthermore, it was found that when V2O3 
is epitaxially grown on LiTaOs, the transition tempera- 
ture is enhanced from the bulk value by 20 K.— Given 
the possibilities of compressive or tensile forces acting on 
the V2O3 due to the epitaxial relationship of the V2O3 
impurity with the CaV204 host and the possibility of 
nonstoichiometry of the V2O3 impurity phase, one can 
see how the transition temperature of the V2O3 might be 
depressed or enhanced from the bulk value by « 30 K as 



we found for the coherently grown V2O3 impurity phases 
in our two annealed crystals in Sec. Ill CI 

In summary, then, it appears that there is an intrinsic 
phase transition in the two annealed CaV204 crystals at 
about 200 K that has no obvious source. We speculate 
that this transition may be the long-sought chiral phase 
transition originally postulated by Villain,^ where there 
is long-range chiral order but no long-range spin order be- 
low the transition temperature, and the long-range chiral 
order is lost above the transition temperature. 

B. Magnetic Susceptibility and Magnetic Heat 
Capacity 

In separate experiments to be described elsewhere;^ 
we have carried out inelastic neutron scattering mea- 
surements of the magnetic excitation dispersion relations 
for CaV204 single crystals. We find that the dispersion 
along the c-axis (in the vanadium chain direction) is sig- 
nificantly larger than in the two perpendicular directions. 
Above the Neel temperature Tn, the magnetic suscepti- 
bility in Fig. [TUj is nearly isotropic. Thus a quasi-one- 
dimensional Heisenberg model appears to be appropriate 
for the spin interactions in CaV204. 

The crystal structure suggests the presence of spin 
5 = 1 zigzag spin chains along the orthorhombic c-axis. 
We report here exact diagonalization (ED) calculations 
of the magnetic spin susceptibility versus temperature 
x(r) and the magnetic heat capacity C(T) of spin 5 = 1 
J1-J2 Heisenberg chains containing = 8, 10, and 12 
spins for J2/J1 ratios from —1 to 5, and containing 14 
spins for J2/J1 = 0. We also report the results of quan- 
tum Monte Carlo (QMC) simulations of x(T) and C{T). 
These simulations were carried out with the ALPS di- 
rected loop application'^- in the stochastic series expan- 
sion framework'^^ for chains with A'^ = 30 and 60 spins 
and J2/J1 = 0. Here Ji and J2 are the nearest- neighbor 
and next-nearest-neighbor interactions on a linear chain, 
respectively. The spin Hamiltonian is the A = 1 special 
case of Eq. ([1]), given by 

N 

"H = ^( JiSj • Si+i + J2S1 ■ Si+2), (3) 

i=l 

where S is a spin-1 operator. Periodic boundary condi- 
tions are imposed, so the chains become rings. Ji is al- 
ways positive (antiferromagnetic) whereas J2 was taken 
to be either positive or negative (ferromagnetic). This 
chain is topologically the same as a zigzag chain in which 
Ji is the nearest-neighbor interaction between the two 
legs of the zigzag chain and J2 is the nearest-neighbor in- 
teraction along either leg of the zigzag chain. For J2 = 0, 
the A^ spins are all part of the same nearest-neighbor ex- 
change (Ji) chain. For Ji = 0, two independent isolated 
equivalent chains are formed, each containing A^/2 spins 
and with nearest-neighbor exchange J2. This effect can 
be quantified using the T = correlation length ^ which 
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FIG. 17: Calculated magnetic spin susceptibility x for spin 
S — 1 J1-J2 Heisenberg chains containing A'' = 12 spins versus 
temperature T, where Ji and J2 are the nearest-neighbor and 
next-nearest-neighbor exchange interactions in the chain. The 
curves from top to bottom on the right are for J2/J1 = 0, 0.2, 
0.4, 2.0, 2.5, 5.0. 



FIG. 18: Magnetic susceptibility x versus temperature T for 
the spin S = 1 J1-J2 Heisenberg chain containing A'' = 12 
spins. Here, Jmax = max(Ji, J2). Pairs of curves for J2/J1 
ratios that are reciprocals of each other become the same at 
high temperatures. 



has been computed in Ref. |37|. We find that we can reach 
ratios N/^ which are at least 2 for J2/J1 ^ 0.6 whereas 
^ becomes comparable to or even bigger than the sys- 
tem sizes N which are accessible by ED for larger J2 / Ji . 
Accordingly, our finite chain data become a poorer ap- 
proximation to the infinite J1-J2 chain for large J2/J1. 
This is exemplified below in Fig. [20] where the data for 
chains containing different numbers N of spins exhibit 
an increasing divergence from each other with increasing 
J2/J1. 

We will compare the spin susceptibility calculations 
with the experimental magnetic susceptibility data to es- 
timate the Ji and J2/J1 values in the J1-J2 chain model 
for CaV204- These values will also be used as input 
to compare the calculated magnetic heat capacity versus 
temperature with the experimental heat capacity data. 



1. Magnetic Susceptibility 

The calculated magnetic spin susceptibility x(T) data 
for the spin S — I J1-J2 Heisenberg chain model are in 
the dimcnsionless form 



xJi 



kBT 
Ji ' 



(4) 



where N is the number of spins, g is the spectroscopic 
splitting factor (5-factor) of the magnetic moments for 
a particular direction of the applied magnetic field with 
respect to the crystal axes, fiB is the Bohr magneton, 
and /cb is Boltzmann's constant. Calculated x(^) data 
sets for = 12 were obtained by exact diagonalization 
assuming periodic boundary conditions (ring geometry) 
for J2/J1 ratios of —1, —0.8, 2.0, 2.5, 5. Exam- 
ples of the calculations for a selection of J2 / Ji values are 



shown in Fig. 1171 Each chain has an energy gap (spin 
gap) from the nonmagnetic singlet ground state to the 
lowest magnetic excited states. '^^ No interchain (between 
adjacent zigzag chains) interactions are included in the 
calculations. These calculations are not expected to ap- 
ply to our system at low temperatures where we see long- 
range antiferromagnetic ordering. However, we expect to 
be able to obtain approximate estimates of Ji and J2 by 
fitting the observed susceptibility data around the broad 
peak in the susceptibility at ~ 300 K. 

At high temperatures ^ Jmax, where Jmax = 

max(Ji, J2), one expects that the calculated xJmax ver- 
sus fce^^/Jmax should be nearly the same upon inter- 
change of Ji and J2, i.e., the same for pairs of J2/J1 
ratios that are reciprocals of each other. [This is be- 
cause all spins in the zigzag chain are equivalent, and 
at high temperatures the Curie- Weiss law is obtained. 
The Weiss temperature 6 only depends on the numbers 
of nearest neighbors 2; to a given spin and the corre- 
sponding interaction strengths J {0 ^ ziJi + Z2J2 with 
zi = Z2 = 2), which is invariant upon interchange of Ji 
and J2 .] This expectation is confirmed in Fig. [TH] where 
such plots are shown for J2/J1 = 1/5 and 5; 1/2 and 
2; and 1. The data for J2/J1 = 1/2 and 2, and for 
J2/J1 — 1/5 and 5, are seen to be about the same for 
temperatures /cbT/ Jmax ^ 4, respectively. 

The experimental magnetic susceptibility data of 
CaV204 will be fitted below by the calculated suscep- 
tibility xiT) of a single 5" = 1 chain (J2/J1 = 0). Such 
integer-spin chains are known as Haldane chains.'^* We 
will therefore test here the sensitivity of the calculations 
to the number of spins N in the chain for this fixed 
J2/J1 value. Shown in Fig. [19] are exact diagonalization 
(ED) calculations of x(T) for J2/ Ji = and A = 12 
and 14, and quantum Monte Carlo (QMC) simulations 
for J2/J1 = and A^ = 60. On the scale of the fig- 
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FIG. 19: Magnetic susceptibility x versus temperature T cal- 
culations for the spin 5* = 1 Heisenberg chain with nearest 
neighbor exchange interaction Ji and next-nearest-neighbor 
interaction J2 — 0. The calculations were carried out using 
exact diagonalization (ED) for A*' = 12 and 14, and by quan- 
tum Monte Carlo (QMC) for N = 60. 




FIG. 20; The product y^'^^^T"'^^ versus J2/J1 for the spin 
S — 1 J1-J2 Heisenberg chain containing A'' = 8, 10, 12, 30 or 
60 spins, where we have assumed A*' — 2Na and g = 1.97. The 
data for N = 8, 10, and 12 were obtained using exact diago- 
nalization calculations. The data for A'' = 30 and 60 were 
obtained from quantum Monte Carlo simulations; most of 
these data were obtained for A' — 60, except for J2/J1 — 0.2, 
1.8 and 2 where we used A'^ = 30. The horizontal dashed 
line is the experimental value of ^'^^^T'"^^ for CaV2 04 from 
Eq. (I?)). Comparison of this experimental value with the cal- 
culations indicates that within the J1-J2 model, J2/J1 ^ 
(or J1/J2 0) in CaV204. 



ure, the results of the three calculations can hardly be 
distinguished. These data are fully consistent with pre- 
vious transfer-matrix renormalization-group results for 
x(T)^ii£l In Table Ell the values of the maxima in the 
magnetic susceptibility ^'^'^^ and also of the magnetic 
heat capacity C"^'^^ (see below) and the temperatures 
jimax g^^^ jimax which they respectively occur are 
listed for the different calculations. For all three cal- 



culations, the maximum in the susceptibility occurs at 
about the same temperature ksT^^^/Ji « 1.30, which 
may be compared with previous values of 1.35 (Ref. l4ll ) 
and 1.32(3).^^ Probably the most accurate values for the 
susceptibility are those of Ref. 0, as listed in Table IVII 
From the theoretical x(T) data, for each value of J2/ Ji 
one can obtain the value of the normalized temperature 
k-gT™^^ / Ji at which the maximum in the susceptibil- 
ity occurs, and the normalized value of the susceptibility 
^max j^^^^2^| at the maximum. For a given value of 
J2I the product of these two values is a particular 
dimensionless number 



X 



maxT^max 



(5) 



that does not contain either exchange constant. 

The spectroscopic splitting tensor (g-tensor) for vana- 
dium cations is found to not depend much on either the 
oxidation (spin) state of the V cation or on its detailed 
environment in insulating hosts. The physical origin of 
this insensitivity is the small magnitude of the spin-orbit 
coupling constant for the vanadium atom. Typical val- 
ues for the spherically-averaged g-factor (g) are between 
approximately 1.93 and 1.97, with the individual com- 
ponents of the diagonal g-tensor lying between 1.90 and 
2.00. For example, for V"*"^ in single crystals of AgCl, 
one obtains (g) = 1.970(3);^^ for V"*"'^ in guanidinium 
vanadium sulfate hexahydrate, (5) = 1.94(1))^'' for V"'"'' 
in Ti02, (.g) = 1.973(4)4^ 

On the basis of the above discussion we set g ~ 1.97 for 
the spin S" = 1 in Eq. (O. Then setting N = 2N p,, 
where TVa is Avogadro's number and the factor of 2 comes 
from two atoms of V per formula unit of CaV204, the 
expression in Eq. ([5]) becomes 



X 



max 2^ in ax 



2.91 cm3K/mol 



(6) 



where a "mol" refers to a mole of CaV204 formula units. 
Then from Eq. ([6]) and the calculated x(T) data for dif- 
ferent values of J2/J1, the calculated -^^^^T"^"^^ versus 
J2/J1 for CaV204 is shown in Fig. [201 From the figure, 
the dependence of T^maxymax J^/Ji is about the same 
for N = 8, 10, and 12 for J2/J1 < 0.6 which is consistent 
with a short correlation length ^ ^ 6 for < J2/«/i < 0.6 
(see Ref. IstI ). However, the curves for the different val- 
ues of N are quite different at larger values of J2/J1; 
the behavior versus N even becomes nonmonotonic in 
this parameter region. The QMC results also shown in 
Fig. [20] nevertheless indicate that the ED calculations for 

= 12 sites yield a good approximation to the infinite N 
limit also for J2/J1 > 1-8. Unfortunately, the QMC sign 
problems are so severe in the region 0.2 < J2/J1 < 1.8 
that here we cannot resolve the maximum of x with our 
QMC simulations. 

The experimental susceptibility Xcxp 

(T) data for 

CaV204 in Fig.lTUlcan be written as the sum Xcxp{T) = 
x(T) +X01 where x(r) is the spin susceptibility (which is 
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TABLE VI: Calculated values of the maxima in the mag- 
netic spin susceptibility x™^'' ^^nd magnetic heat capacity 
C"""' and temperatures TJ"'''' and T^'''' at which they oc- 
cur, respectively, for the linear spin S = 1 Heisenberg chain 
(Haldane chain) with nearest-neighbor exchange interaction 
Ji and next-nearest-neighbor interaction J2 = 0. The re- 
sults of exact diagonalization (ED) and quantum Monte Carlo 
(QMC) calculations are shown. Here A'' is the number of spins 
in the chain, g is the g-factor, /xb is the Bohr magneton, and 
is Boltzmann's constant. Also included are the results of 
Ref.l40l. which are probably the most accurate values currently 
available for the susceptibility. 
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FIG. 21: Calculated magnetic susceptibility x versus temper- 
ature T for the J1-J2 chain with Ji/ks ~ 230 K, J2 =0, and 
g = 1.97, and with a temperature- independent orbital contri- 
bution xo ~ 0.4 X 10"'^ cm'^/mol (solid curve). The calculated 
spin susceptibility at T = is zero, so the zero-temperature 
value of the calculated solid curve is xo- Experimental data 
for annealed CaV2 04 crystal 2-50-cl from Fig. [10] are also 
shown. Comparison of these data for T — > with the calcu- 
lated curve shows that the spin susceptibility along the easy 
&-axis of CaV2 04 is rather large for T — > 0. 

the part calculated above) and xo is the temperature- 
independent orbital susceptibility. From the data in 
Ref. m for V2O3, we estimate xo ~ 0.4 x 10"'^ cm'^/mol 
for CaV204. From Fig. [TO] we then obtain the experi- 
mental value for the spin susceptibility at the maximum 
^max ^ 2.2 X 10"'^ cm'^/mol and for the temperature at 
the maximum T^^"" « 300 K, yielding for CaV204 

^maxymax _ q gg ^^^^3 J^/^^l (7) 

Comparison of this value with the theoretical spin suscep- 
tibility data in Fig. [2D] yields J2/J1 « (or J1/J2 « 0). 
This ratio of J2 / Ji is quite different from the value of 



unity that we and others"'^'''^^ initially expected. The tem- 
perature T^'^'^ « 300 K, combined with k^T^'''' / Ji « 
1.30 from Table [VTl yields Ji/fce = 230 K. Although the 
numerical results shown in Fig. [20] are the least accurate 
in the vicinity of J2/J1 « 1, it seems rather unlikely that 
the value of T^maxjimax oijtained from the J1-J2 chain 
model in the region 0.6 < J2I J\ ^1-8 could be consis- 
tent with the value in Eq. ([7]) expected for CaV204. 

The calculated total susceptibility versus temperature 
for Ji/fcs = 230 K, J2 = and xo = 0.4 x lO'^ cmVmol 
is shown in Fig. [2U Also shown are the experimental 
susceptibility data for annealed CaV204 crystal an-2-50- 
cl from Fig. [TO] where an excellent fit of the average 
anisotropic x(^) data near 300 K is seen. 

Above the Neel temperature, one sees from Figs. [TO] 
and [21] that the susceptibility is nearly isotropic. The 
relatively small anisotropy observed can come from 
anisotropy in the orbital Van Vleck paramagnetic sus- 
ceptibility, from (7-anisotropy arising from spin-orbit in- 
teractions, from single-ion anisotropy of the form DS\ 
E{Sl — Sy), and/or from anisotropy in the spin exchange 
part of the spin Hamiltonian. The relative importances of 
these sources to the observed susceptibility anisotropics 
are not yet clear. The experimental data below 200 K 
in Fig. [21] increasingly deviate from the fit with de- 
creasing temperature. This suggests that other interac- 
tions besides Ji and J2 and/or the presence of magnetic 
anisotropics may be important to determining the spin 
susceptibility above Tn in CaV204. 

For coUinear antiferromagnetic (AF) ordering, one 
nominally expects the spin susceptibility along the easy 
axis to go to zero as T — > 0. Comparison of the the- 
oretical curve with the experimental easy axis (5-axis) 
data Xb{T) in Fig. [^indicates that the zero-temperature 
6-axis spin susceptibility is not zero, but is instead a 
rather large value xT'"\T ^ 0) w 0.9 x 10"^ cm^/mol. 
This finite spin susceptibility indicates either that the 
spin structure in the AF state is not coUinear, that not 
all vanadium spins become part of the ordered mag- 
netic structure below Tn, and/or that quantum fluctu- 
ations are present that induce a nonzero spin susceptibil- 
ity. Such quantum fluctuations can arise from the low- 
dimensionality of the spin lattice and/or from frustration 
effects. As discussed in the Introduction, our recent NMR 
and magnetic neutron diffraction experiments on single 
crystal CaV204 indicated that the magnetic structure at 
4 K is noncoUinearfiii^ which can at least partially ex- 
plain the nonzero spin susceptibility along the (average) 
easy 6-axis at low temperatures. In addition, the reduc- 
tion in the local ordered moment 1.0-1.6 /j,b/(V atom) of 
the ordered vanadium spins found in these studies from 
the expected value gSfiB = 2 /xb/ (V atom) suggests that 
quantum zero-point spin fluctuations could be strong and 
could contribute to the large finite Xb^'" (^^ ^ 0) . 
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FIG. 22: Magnetic heat capacity C versus temperature T for 
the spin S — 1 J1-J2 Heisenberg chain with J2/J1 values from 
— 1 to 5, calculated using exact diagonalization with — 12. 
Here A'' is the number of spins, /cb is Boltzmann's constant, 
and Ji > and J2 are the nearest-neighbor and next-nearest- 
neighbor exchange interactions, respectively. 
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FIG. 23: Magnetic heat capacity C versus temperature T for 
the S = 1 J1-J2 chain with J2/J1 ~ ("Haldane chain"), 
calculated using exact diagonalization (ED) with = 10 and 
12, and quantum Monte Carlo (QMC) simulations for — 60. 



2. Magnetic Heat Capacity 

The magnetic heat capacity C versus temperature T 
was calculated by exact diagonalization for = 12 spins 
S = 1 over the range —1 < J2/J1 < 5. Represen- 
tative results are plotted in Fig. [^H The variation in 
C{T) with N is illustrated in Fig. [23] for J2/J1 = and 

= 12 and 14 from exact diagonalization calculations, 
and for A'' = 60 from quantum Monte Carlo simula- 
tions. The data for the different N are seen to be nearly 
the same. The values of the maxima C'^^^ in the mag- 
netic heat capacity and the temperatures T^'^^ at which 
they occur are listed above in Table IVTl Our results for 
the specific heat are consistent with previous transfer- 
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FIG. 24: Magnetic heat capacity C versus temperature T for 
the spin S = 1 J1-J2 chain with J2/J1 = and Ji/ks = 
230 K, calculated using exact diagonalization with A'^ = 14 
(solid curve) from Fig. 1231 The data points are the AC{T) 
data for annealed crystal an-2-50-cl from Fig. 1151 



matrix renormalization-group computations.— The 
two transfer-matrix renormalization-group results differ 
at high temperatures. Our QMC results for C obtained 
from rings with = 60 sites are in better agreement with 
the older results which apply to the infinite A^ limili^ 
than the more recent results obtained for open chains 
with A^ = 64 sites4a The C(T) data in Ref. \M were 
calculated from a numerical derivative which resulted in 
systematic errors in the data at high temperatures 4i 

We cannot confidently derive the exchange constants 
in CaV204 from fits of our heat capacity data by the 
theory. Extraction of the magnetic part of the experi- 
mental heat capacity is tenuous because of the presence 
of the orthorhombic to monoclinic structural transition 
at Ts « 150 K and the transition(s) at Tsi « 200 K. 
Therefore the relationship of the heat capacities AC (T) 
in Figs. [Ufa) and [15] to the magnetic heat capacities of 
the samples is unclear. Here, we will just compare the 
theoretical magnetic heat capacity calculated for the ex- 
change constants Ji = 230 K and J2 = 0, that were 
already deduced in the previous section, with the exper- 
imental AC(T) data to see if theory and experiment are 
at least roughly in agreement. This comparison is shown 
in Fig. [M] for annealed CaV204 crystal an-2-50-cl from 
Fig. [15] Overall, the theory and experimental data have 
roughly the same magnitude, but the data are system- 
atically below the theoretical prediction. This is likely 
caused by the heat capacity of the nonmagnetic reference 
compound CaSc204 being somewhat different from the 
lattice heat capacity of CaV204. We note from Fig. [13] 
that a difference of 5 J/mol K between the theory and 
experiment in Fig. [M]is only about 5% of the total heat 
capacity of the samples at 200 K. In addition, we have 
not included in the theory interchain couplings that lead 
to long-range antiferromagnetic order, or the effect of the 
magnetic ordering on the heat capacity including the ef- 
fect of the energy gap in the spin wave spectrum below 
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C. Interchain Coupling 

Within the S = I J1-J2 Heisenberg spin chain model, 
we found above that J2/J1 ~ and Ji « 230 K in 
CaV2 04 near room temperature. Thus the crystallo- 
graphic zigzag vanadium chains in CaV204 act hke 5 = 1 
hnear spin chains with nearest-neighbor interaction Ji. 
This is a so-called Haldane chain"^^. with a nonmagnetic 
singlet ground state and an energy gap for spin excita- 
tions given b}42ii^ ss 0.4105Ji. An interchain coupling 
J± must be present in order to overcome this spin gap 
and induce long-range antiferromagnetic ordering at Tn- 
Pedrini et al^^^ have recently estimated the dependence 
of Tn/Ji on Jj_/Ji using a random-phase approximation 
for the interchain coupling for 5 = 1 Haldane chains. Us- 
ing our values Tn = 68 K, Ji = 230 K and T'"^'^ = 300 K, 
we obtain J±/Ji w 0.04 and J± — 5-10 K. However, 
it should be emphasized that the treatment of Refs. '43 
and [13 assumes a nonfrustrated interchain coupling ge- 
ometry such that the result Jx « 10 K should be con- 
sidered as a lower bound. Still, the value of J±/Ji is 
sufhciently small that a redetermination of J2/J1 and 
Ji, from a Ji-J2-J± model fitted to the observed suscep- 
tibility data for CaV2 04 near room temperature, would 
yield very similar values of J2/ Ji and Ji to those we have 
already estimated using the isolated chain J1-J2 model. 
Additional and more conclusive information about the in- 
terchain coupling strength(s) will become available from 
analysis of inelastic neutron scattering measurements of 
the magnetic excitation dispersion relationsJ^ 



V. SUMMARY 

We have synthesized the 5* = 1 spin chain com- 
pound CaV2 04 in high purity polycrystalline form and 
as single crystals. Our magnetic susceptibility x{T) and 
ac magnetic susceptibility Xac(r) measurements do not 
show any signature of a spin glass-like transition around 
20 K that was previously reportedjiSiii We instead ob- 
serve long-range antiferromagnetic ordering at sample- 
dependent Neel temperatures Tn ~ 50-70 K as shown in 
Table IIVI The Neel temperature and the orthorhombic- 
to-monoclinic structural transition temperature Tg in Ta- 
ble IIVI both show a large systematic variation between 
different samples. Those temperatures for an unannealed 
crystal are each less than those for an annealed crystal 
which in turn are less than those for a sintered polycrys- 
talline sample. The cause of these large temperature dif- 
ferences, especially between as-grown and annealed single 
crystals, is unclear. The transition temperature differ- 
ences may arise from small changes in oxygen stoichiom- 
etry (< 1 at.%, below the threshold of detection by TGA 
or XRD) and/or from structural strain, both of which 
may be reduced upon annealing the as-grown crystals at 



1200 °C in 5% H2/He. In addition, other small chemical 
differences and/or structural defects may be relevant. 

Our heat capacity Cp(T), linear thermal expansion 
a(T), and x(^) measurements reveal distinct features 
at the orthorhombic-to-monoclinic structural transition 
temperature Ts identified from our diffraction studieSfi^ 
We inferred from a combination of structural studies and 
physical property measurements that the origin of the 
third transition at Tsi ~ 200 K in one of our annealed 
crystals was mostly due to the structural transition in 
the V2O3 impurity phase that grew coherently upon an- 
nealing the crystal. In the other annealed crystal, we 
ruled out this source and we are thus left with a transi- 
tion at Tgi with unknown origin. We speculate that this 
transition may be the long-sought chiral phase transition 
originally postulated by Villain in 1977.^ 

The x(^) shows a broad maximum at about 300 K 
indicating short-range antiferromagnetic (AF) order- 
ing in a low-dimensional antiferromagnet as previously 
observediSiii and the x{T) above Tn in single crystals is 
nearly isotropic. The anisotropic x{T) below Tn shows 
that the (average) easy axis of the antiferromagnetic 
structure is the orthorhombic &-axis. The magnetic spin 
susceptibility along this axis is found to be a large finite 
value for T — > 0, instead of being zero as expected for 
a classical collinear antiferromagnet. This result is con- 
sistent with our observed noncollinear magnetic struc- 
ture below Tnj ^^'^^ In view of the fact that CaV204 is a 
low-dimensional spin system, quantum fluctuations could 
also contribute to both the observed reduced zero tem- 
perature ordered moment and the relatively large zero 
temperature spin susceptibility. 

We analyzed the xC^) data near room temperature in 
terms of theory for the 5=1 J1-J2 linear Heisenberg 
chain, where Ji(J2) is the (next-)nearest-neighbor inter- 
action along the chain. We obtain Ji/fce ~ 230 K, but 
surprisingly J2/J1 ~ (or J1/J2 ~ 0), so the exchange 
connectivity of the spin lattice appears to correspond to 
linear 5 = 1 Haldane chains instead of zigzag spin chains 
as expected from the crystal structure. This result is 
consistent with analysis of our high temperature (up to 
1000 K) magnetic susceptibility measurements on single 
crystal CaV2 04ii^ We estimated here the coupling J± 
between these chains that leads to long-range AF order 
at Tn to be J±/Ji > 0.04, i.e., only slightly larger than 
the value J±/Ji « 0.02 needed^i^ to eliminate the en- 
ergy gap (Haldane gap) for magnetic excitations. 

From our Cp(T) measurements, the estimated molar 
magnetic entropy at Tn is only w 8% of its maximum 
value 2Rln(25 + 1) = 2Rln(3), where R is the molar gas 
constant, and the heat capacity jump at Tn is only a few 
percent of the value expected in mean field theory for 
5 = 1. Both results indicate strong short range antifer- 
romagnetic order above Tn and large values Ji and/or 
J2 > 100 K, consistent with the xC^) data. We also 
compared the Cp(T) data with the theoretical prediction 
for the magnetic heat capacity using the exchange con- 
stants found from the magnetic susceptibility analysis. 
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and rough agreement was found. However, this com- 
parison is not very precise or useful because the struc- 
tural transition at Ts ^ 150 K and the transition(s) at 
Tgi ~ 200 K for our two annealed single crystals, make 
large contributions to Cp(T). In addition, the accuracy 
of the measured heat capacity of the nonmagnetic refer- 
ence compound CaSc2 04 in representing the lattice heat 
capacity of CaV2 04 is unknown. Thus extracting the 
magnetic part of the heat capacity at high temperatures 
from the observed Cp(T) data for comparison with the- 
ory is ambiguous. 

In closing, we note the following additional issues that 
could usefully be addressed in future work. Our analyses 
of our x(^) data for CaV204 to obtain the exchange con- 
stants Ji and J2 were based on fitting the experimental 
x(T) data only near room temperature, since our calcu- 
lations of x{T) all showed nonmagnetic singlet ground 
states, contrary to observation, and could not reproduce 
the observed antiferromagnetic ordering at low temper- 
atures. Calculations containing additional interactions 
(see also below) and/or anisotropics are needed for com- 
parison with the lower temperature data. 

The orthorhombic crystal structure of CaV204 at 
room temperature contains two crystallographically in- 
equivalent but similar V^'^ 5=1 zigzag chains. These 
chains may therefore have different exchange constants 
associated with each of them. For simplicity, our x{T) 
data were analyzed assuming a single type of zigzag 
chain. Furthermore, the extent to which the transitions 
at Tsi and Tg affect the magnetic interactions is not yet 
clear. 

From crystal structure considerations, one expects that 
J2/J1 « 1 in CaV204,-i^'ii instead of J2/J1 « as 
found here. This suggests that additional magnetic inter- 
actions and/or anisotropy terms beyond the Heisenberg 
interactions Ji and J2 and interchain coupling J± con- 
sidered here may be important. In addition to single 
ion anisotropy and other types of anisotropy, we mention 
as possibilities the Dzyaloshinskii-Moriya interaction, bi- 
quadratic exchange, and cyclic exchange interactions 
within the zigzag chains. When such additional terms 
are included in the analysis, the fitted value of J2/J1 
could turn out to be closer to unity. A four-spin cyclic 
exchange interaction has been found to be important to 
the magnetic susceptibility in cuprate spin ladders.^ In 
these spin ladders, there are exchange interactions J and 
J' between nearest-neighbor Cu+^ spins 1/2 along the 
legs and across the rungs of the spin ladder, respectively. 
For the 5 = 1/2 two- leg ladder compound SrCu203, if 
only J and J' are included in fits to the data, one obtains 
J'/ J « 0.5.^ However, by also including the theoretically 
derived cyclic four-spin exchange interaction, the ratio 
J' / J increases from 0.5 to a value closer to unity, as ex- 
pected from the crystal structure. 



Pieper et al. have recently proposed a very different 
and very interesting model to explain the inference that 
J2 / Ji « around room temperature which involves par- 
tial orbital ordering of the two d-electrons of V among 
the three t2g orbitals.^^ Furthermore, in order to explain 
the magnetic structure at low temperatures, they deduce 
that the nature of the orbital ordering changes below Ts 
such that the effective spin lattice becomes a spin-1 two- 
leg ladder. 

It has been well documented that fits of magnetic sus- 
ceptibility data by theory tests only the consistency of a 
spin model with the data, and not the uniqueness of the 
model. A good example of this fact arose in the study 
of the antiferromagnetic alternating exchange chain com- 
pound vanadyl pryophosphate, (VO)2P202, the history 
of which is described in detail in the introduction of 
Ref. [51I. The ultimate arbiter of the validity of a spin 
model is inelastic neutron scattering measurements of the 
magnetic excitation dispersion relations in single crystals. 
Theoretical calculations of the exchange interactions are 
much needed and would also be valuable in this regard. 

Finally, the origin of the intrinsic heat capacity anoma- 
lies at Tsi « 200 K for the two annealed single crystals 
of CaV2 04 needs to be further studied. We speculate 
that this transition may be the long-sought chiral phase 
transition originally postulated by Villain in 1977i^ 

Note added — After this work and this paper were 
nearly completed, Sakurai reported a very interesting 
and detailed study of the magnetic and electronic phase 
diagram of polycrystalline samples of the solid solution 
Cai_a;Naa;V204 prepared under high pressurci^ 
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